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Two techniques are described that simplify the experimental requirements for measuring and 
manipulating quantum information stored in trapped ions. The first is a new technique using 
electron shelving to measure the populations of the Zeeman sublevels of the ground state, in an ion 
for which no cycling transition exists from any of these sublevels. The second technique is laser 
coohng to the vibrational ground state, without the need for a trap operating in the Lamb-Dicke 
hmit. This requires sideband cooling in a sub-recoil regime. We present a thorough analysis of 
sideband cooling on one or a pair of sidebands simultaneously. 
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Laser cooling and electron shelving detection techniques |^-^ were developed in the pursuit of better control 
of a basic physical system, and high precision experiments. Up to now such experiments have concentrated on atomic 
transitions that offer good prospects as frequency standards, or that allow sensitive detection of some physical effect. 
The relatively new subject of quantum information theory has led to interest in a different experimental approach. 
Here, an important aim is to realise the preparation, coherent control, and detailed measurement of some sufficiently 
complex physical system. This allows one to realise simple networks of elementary quantum operations (quantum 
gates), and to test basic ideas in quantum information theory, such as data compression, quantum error correction, 
and simple quantum algorithms. The new feature is that such experiments do not require high precision as an end in 
itself (unlike frequency standards), although high precision is of course very desirable, nor do they require sensitivity 
to some physical effect. Instead, one's attention is focused purely on the logical properties of the state transformations 
and measurements that can be carried out. 

The significance of experiments on trapped ions to quantum information physics was recognised by Cirac and ZoUer 
[0, who proposed the use of a linear ion trap to realise the essential elements of a quantum computer. Although a 
real ion trap will not permit large-scale quantum computing [pl]-p^ , such a system should allow experiments on a few 
tens of qubits, and is one of the most promising for such a purpose. Currently a few groups worldwide are developing 
experiments to pursue such ideas fl^ . 

In this paper we consider the question "what is the simplest method to build and operate a linear ion trap quantum 
information processor?" One of the available options is to use an ion such as calcium or strontium whose ground state 
has total spin half. The use of the two Zeeman sublevels of the ground state of each ion naturally suggests itself as 
a means to store each qubit, but currently this option has not received much attention, partly because the method 
to measure the final state, the "readout" in computer terminology, is not obvious. We propose in sections 1 and 2 
a simple electron shelving technique that will permit such measurements. Our method does not realise an ideal von 
Neumann measurement (projection onto an orthonormal basis) but nevertheless permits efficient quantum computing 
and quantum state tomography. This is discussed in section 3. 

In section 4 we consider state preparation in this system. The most difficult part of the state preparation is the 
cooling of the motional degrees of freedom. We present the results of a study of sideband cooling in the regime where 
the confinement is of intermediate tightness (Lamb-Dicke parameter of order one) , since this is typically the regime 
in which one would wish to operate the trap as a "processor" . We find that first sideband cooling operates well for 
Lamb-Dicke parameters up to 77 ~ 0.6, and indeed sub-recoil cooling is possible. By combining excitation of the first 
sideband with a higher one, cooling to the ground state can be achieved for Lamb-Dicke parameters up to 77 ~ 3. We 
find that the best choice of sideband is higher than one might expect, and the analysis reveals the optimum ratio of 
laser intensities on the two sidebands. 



I. CHOICE OF TRANSITION 



Consider a linear ion trap, in which the vibrational frequencies along the three axes, for a single confined ion, are 
LOx — ujy ^ oJz- Each ion must store one qubit (or possibly more than one), so we require two long-lived states in 
the ion, between which transitions can be driven coherently. The recoil energy is defined Ej^ = {fik)"^ / {2M) where k 
is the wavevector of the radiation used to drive the transition and M is the mass of a single ion. The Lamb-Dicke 
parameter is 

ri = mcos{e){ER/huJz)^''^ (1) 
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where 9 is the angle between the laser beam(s) and the z axis (k-z = k cos{0)) and m = 1 for single photon transitions, 
m = 2 for Raman transitions, assuming the geometry ki • z = —\s.2 ■ z. 

We require the Lamb Dicke parameter not to be vanishingly small, in order that the radiation can affect the motional 
state of the ions, and therefore optical rather than radio-frequency radiation must be used. Since the separation of 
adjacent trapped ions is in the region 10 to 100 /im, radiation of sub-micron wavelength is also preferable in order to 
address individual ions. These considerations lead to two approaches to information processing experiments: either a 
single-photon transition between metastable states separated by optical frequencies is adopted, or an rf transition is 
driven using optical radiation by means of the stimulated Raman effect. 

We will consider the case of Raman transitions, since the laser frequency stability requirement is much less restrictive 
[0 . The Raman transition is sensitive only to the difference of the frequency of the two laser beams used to drive 
it. If both beams are derived from the same laser, the laser's frequency fluctuations cancel out, and we only have to 
concern ourselves (to first order) with power fluctuations. A single-photon transition, by contrast, would require a 
laser system close to the state of the art in frequency stability. 

Atomic energy levels that have long natural lifetimes and separations in the rf regime almost always owe their 
energy separation to the hyperfine interaction, or the Zeeman or Stark effects. The singly-charged ^Be ion is a good 
candidate for information processing purposes because it has hyperfine structure. However, the transition wavelengths 
in ions such as *'^Ca and *^Sr are more easily accessible with current laser technology. This leads us to consider these 
ions as candidates for quantum information experiments. Since they lack hyperfine structure, and we wish to use 
long-lived states with rf separation, we store each qubit in the two-dimensional Hilbert space spanned by the two 
Zeeman sublevels \M = ±1/2)^ of the ground state. The original proposal of Cirac and Zoller |l^ made use of a third 
long-lived state in each ion in order to bring about quantum gates such as XOR, which is not available here. However, 
we can avoid the need for such a state either by use of the "magic Lamb-Dicke parameter" method of Monroe et al. 
I psf , or by using the first excited state of the second normal mode of oscillation. 

One disadvantage of our choice is that the transition \M = —1/2)^ ^ \M = -t-1/2)^ is sensitive to the local magnetic 
field to first order. Although this would rule out such a transition as a frequency standard, it does not for information 
processing. The reason is that we only require short term stability, not absolute precision, in the energy level 
separations. The actual value of these separations is totally immaterial, as long as the Raman resonance condition 
can be found precisely in any given experiment. An advantage of our choice is that there arc no further sublevels in 
the ground state into which population can "leak" , which is one source of error in quantum information processing. 

A major apparent problem with our choice is in the final state detection. In quantum processing, we require a 
sensitive measurement of the final state of each individual qubit. Sufficient signal to noise ratio to make such a 
measurement is quite rare in physics, but it can be done for individual trapped atoms by exciting fiuorescence on 
an optical transition that is resonant for one state (say |A/ = +1/2)^) but not for the orthogonal state. Typically 
at least a few thousand cycles of excitation followed by spontaneous decay must be completed in order to get an 
unambiguous signal. This is possible if the optically excited state can only decay to the \M = -1-1/2)^ ground state, 
that is not to |M = —1/2)^ and not to some other metastable state. However, in ions such as Ca, Sr, Ba and Hg all 
the short-lived excited states (one of which must be used for efficient fluorescence detection) have a chance of decaying 
to a metastable state (see fig. 1), causing fiuorescence to stop after about 3 to 10 cycles. This is the problem we wish 
to address. 



II. ELECTRON SHELVING ZEEMAN STATE MEASUREMENT 

Our scheme to detect the final Zeeman ground state of a single ion of Ca, Sr, or a similar candidate, requires three 
laser wavelengths, exciting the transitions S1/2 ^ P3/2, S1/2 <-> P1/2, and -D3/2 <^ P1/2 (see fig. 1). Since the latter 
two are required for laser cooling in any case, the main cost of the scheme is one additional laser wavelength. The 
wavelengths for '*°Ca are 394 nm, 397 nm and 866 nm respectively. We will use these values as a shorthand to refer 
to the transitions, although obviously the method works essentially unchanged for all similar ions. 

A magnetic field is imposed in a direction along the wavevector of the the 394 nm radiation, which typically is not 
the z axis of the ion trap. We adopt the magnetic field direction as the quantization axis for the purpose of defining 
spin states and polarizations. The 394 nm radiation is polarized tr"*". The polarization of the other beams is almost 
unconstrained, linear polarization at right angles to the magnetic field direction (i.e. equal amounts of cr~^,a~) will 
do. 

During the first part of the detection scheme, the magnetic field must be large enough to enable the Zeeman 
components of the S1/2 P3/2 transition (394 nm) to be resolved. In other words, their splitting is of the order of or 
greater than the inverse of the lifetime of the excited state {P3/2)- The 394 nm radiation is tuned to resonance with 
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the 5*1/2, M = +1/2 ^ P2,/2tM = 3/2 transition, and therefore excites transitions from \A4 = +1/2)^ much more 
strongly than from \M — — 1/2)^. 

Let n± be the initial populations of the ground state sublevels, i.e. the diagonal elements of the single-ion density 
matrix in the \M — ±1/2)^ basis. 

The measurement proceeds in two steps: shelving followed by fluorescence. First the 394 nm radiation is pulsed on 
for a time sufficient to cause about 10 cycles of excitation and spontaneous decay. After this time the population that 
was originally in \M = +1/2)^ has been optically pumped (shelved) into the D states, most of it (about 90%) in the 
state, while the population in |M = — 1/2)^ is almost unaffected. For a single ion experiment, this means that 
the state jumps abruptly to the D manifold or to \M = —1/2) with probabilities proportional to 11+, n„ respectively. 
Next, the 394 nm radiation is turned off, the magnetic fielcf is reduced, and the two other wavelengths are turned 
on. Fluorescence is detected during a time equal to a few percent of the lifetime of the I?5/2 state. Since this time 
is more than 10^ times longer than that of the P1/2 state, of the order of a million cycles of fluorescence can take 
place, making a large accumulated signal. The probability of observing fluorescence is proportional approximately to 
n_ + 0.in+. A further signal for normalisation or checking purposes can be obtained after the U5/2 state decays, if 
the ion state jumped there during the shelving pulse. 

We have analysed this scheme by numerically solving the complete set of rate equations, including the 18 Zeeman 
sublevels and 40 allowed electric dipole transitions in the problem. This enables us to calculate the requirements on 
magnetic field strength and laser pulse timing, and to interpret the signal precisely. A rate equation rather than optical 
Bloch equation approach is sufficient since during the shelving pulse, optical coherences play no role in determining 
populations, and during the fluorescence detection we impose a magnetic field sufficient to prevent accumulation of 
population in a "dark state", that is, a spin state in the -D3/2 Zeeman manifold that does not couple to the 866 nm 
radiation. For this, the dark state must Larmor precess to a bright state much faster than the optical pumping rate 
among the Zeeman sublevels of the D3/2 state. The calculation tells us this optical pumping rate. 

In order to write the equations, let us number the five energy levels in order of increasing energy, as in fig. 1. Let 
iV^J^ be the population of the M Zeeman sublevel of the nth level. The total rate of change of population n'^'^ is 

a sum of spontaneous and driven terms. The spontaneous feeding and decay terms are 
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where C'jjn) m ^""^ ^^"^ Clebsch-Gordan coefficients. 

The driven terms are optical excitations caused by the radiation field, taking into account the Rabi frequency fJ, 
detuning 5 and the polarisation of the radiation. For example, during the initial "shelving" stage when only cr+ 
polarized 394 nm radiation is present, tuned to resonance with the Zeeman-shifted Si/2,M — 1/2 ^ P3/2, A/ — 3/2 
transition, the driven terms are 

n%2~n[%) ^l^hg^il^^BB/h), (4) 



dt ^ 

' dr 



dt 



dr 

where the lineshape function is defined by 

T^,=Y.rji^ (6) 
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We obtain the spontaneous decay rates Tji from ||l^. For Ca, they are = 27r x 20.7 MHz, r24 = 27r x 1.69 MHz, 
Fis = 27r X 21.5 MHz, T25 = 27r x 0.177 MHz, Fgs = 27r x 1.58 MHz. The Clebsch-Gordan coefficients and Lande g 
factors are obtained by calculation in the LS coupling approximation. Note that we have neglected the far-off-resonant 
excitation rates in the rate equations written above. They were in fact included in the numerical calculations, though 
their influence is negligible. 
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Wc can quantify the success of the spin measurement by two parameters: the abihty to discriminate the \M = ±1/2)^ 
states, and the size of the signaL We parametrize the former by calculating a "shelving efficiency" e defined to be 
e = S-^- — S-, where 5+ = ^m' population in the -D5/2 levels after the shelving pulse, when the population 

is initially all in \M = +1/2)^ (i.e. 11+ = 1), and S- is the same quantity when the population is initially all in 
\M = —l/2)g (so n_ = 1). In fig. 2 we show e as a function of the magnetic field B, for two values of the Rabi 
frequency, f2i5 = (dashed line) and f2i5 = Ti^/A (full line), with the shelving pulse duration optimized for 
maximum e at each value of B. At large fields, e tends to r35/(r35 + r25) 0.899. The field required to obtain 
e > 0.5 is of order S ~ 3 x 10"'^ T, which is easily produced in the lab. 

It is important to get the shelving pulse duration right, since if it is too long then all the ground state population 
will be pumped to the D states, irrespective of the initial values of II±. This duration, and the required magnetic 
field, can be estimated by a reduced set of rate equations, in the limit that all rates are fast compared to two, namely 
the rate Ri of transfer from \M = -1/2)^ to \M = +1/2)^, and the rate R2 of transfer from \M = +1/2)^ to the D 
states. These are 

Ri = l^l^h95{ll^BB/h), (7) 



The reduced set of equations are the following: 



= -R,N%„ (9) 
= R,N%- R,N^^;i^, (10) 

R2N[% (11) 
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where the final equation accounts for all population in the D levels. Writing 5"+ = 0.899iV(^)(i) when n'^^^}^{0) 
S- = 0.899iV(^)(t) when ivi\^/2(0) = 1> we obtain the solution 



e{t) = 0.899—5^ (e-^i* - e"-^^*) . (12) 
R2 — Ri 

This solution reproduces the result of the complete calculation quite well, see fig. 3. The optimum pulse length, when 
de/dt = 0, is therefore 

_ ln(fi2/fii) 

Imax — jj D ' ^ ) 

-K2 — ill 

and the requirement on the magnetic field is R2 » i?i which is roughly equivalent to the requirement that the Zeeman 

splitting is large compared to Fs, as one would expect. The typical order of magnitude for tmax is lO/Fss ~ 1 /is. 

To complete the analysis, we also solved the complete set of rate equations during the detection phase, when laser 
radiation at 397 nm and 866 nm is present (polarised perpendicular to the magnetic field direction). The observed 
fluorescence rate will be proportional to the population of the P1/2 level, N^^l^2 + ^+1/2 • show in fig. 4 the steady 
state value of this population as a function of magnetic field, with the lasers always tuned to the zero-field resonance 
frequencies, and initial condition n_ = 1. Comparing fig. 4 with fig. 2, the conclusion is that it is not absolutely 
necessary to reduce the magnetic field between the shelving and detection phases, as long as sufficient laser power is 
used during the detection phase, but that a reduction to 10""^ T is useful to increase the observed signal. Since the 
£>5/2 state lifetime is long, there is plenty of time to reduce the magnetic field if so desired. 

We find also that the steady-state P1/2 population remains of order 1/4 when the magnetic field is large enough to 
produce Larmor precession much faster than the optical pumping rate among the Zeeman sublevels of -D3/2- This is 
made possible by the favourable condition r24, and is necessary in order to avoid population accumulating in 

a dark state of D3/2- 
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III. REQUIREMENTS ON MEASUREMENTS FOR QUANTUM COMPUTING 



The scheme we have discussed does not achieve an ideal measurement, which would be represented by the value 
6 = 1. In this section we will discuss whether this is a problem. 

First consider measuring a single ion. Any e > will allow the populations II± to be deduced by repeating the 
measurement many times (preparing the same state each time). It is sufficient that the graph of expected signal 
(the mean number of times fluorescence is seen during the detection phase) verses II-i. is single- valued, since then n_|_ 
(and therefore n_) can be deduced unambiguously from the mean signal. If e is small the measurement will need 
to be repeated more times in order to get good statistics to estimate 11+ , but our value e ~ 0.9 is not significantly 
worse than e = 1 in this respect. Methods available for other ions, such as probing a cycling transition in Be, have e 
closer to 1, but the signal to noise ratio is worse because imperfection in the laser polarisation limits the number of 
fluorescence cycles per measurement to a few thousand, instead of about a million as in our method. 

Now consider a many-ion quantum processor containing N ions. It is clear that e ~ 0.9 is sufficient for information 
processors of small or intermediate size, simply because remains non- negligible (it falls to 0.01 at = 44). The 
question is whether such non-ideal measurements remain useful for a large quantum computer, containing possibly 
thousands of qubits. Although ion traps may not offer the best technology to build large quantum computers in the 
future, the question of measurement errors will remain significant. We will discuss this using the language of ion 
traps; the results can be translated directly to other systems. 

In the context of quantum computation, the final state of the computer represents the output of some algorithm. 
First consider the case that this output (before measurement) is a product state where each ion is either in the state 
\M = —l/2)g or \M = +l/2)g. The probability of obtaining a correct readout after a single run of the algorithm is 

e''^ where A^ is the number of ions in |M = -1-1/2)^. It might be imagined that this exponential dependence on A^ 
renders the complete algorithm (including final measurement) computationally inefficient. In fact this is not the case. 
On repeating the whole algorithm and final measurement r times, independent statistics are gathered on the final 
state of each ion, with the aim of deciding whether the mean number of times fiuorescence is obtained is equal to r 
or r(l — e). The probability that this mean result is interpreted correctly is of order [1 — (1 — e)*"]^. The number of 
repetitions required to make the overall success probability close to one is therefore r ^ 0(log A^/ log(l — e)^^). The 
method is thus efficient. 

In the case that the final state \(j)f) (before measurement) of the quantum computer is not a single product state, 
simply repeating the whole computation may not succeed since the state \4>f) = o-u \u) can consist of a superposi- 
tion of product states \u) in which the number of terms in the superposition is exponentially large, as for example in 
Shor's algorithm. In this case the problem can be avoided by error-correction coding p^ , using a repetition code in 
the measurement basis. The coding requires a repetition r of the same order as that given in the previous paragraph. 
The measurement of the computer's final state must be post-processed (by classical computation) to implement the 
error correction, after which the value u for one of the product states \u) in \(f)f) can be deduced with high probability. 
The situation is now equivalent to having ideal measurements. 

Large quantum computers will rely heavily on quantum error correction, applied repeatedly during the computation. 
Such fault-tolerant methods work best when the error syndromes can be measured, rather than handled by unitary 
quantum networks p^ ], since the non-trivial syndrome interpretation can then be carried out reliably by a classical 
computer. Hence it is important to ask whether e 0.9 renders the syndrome measurement too unreliable. A 
preliminary calculation using the methods of fault-tolerant syndrome extraction suggests that reliable computing can 
be achieved with this value of e, however this is a long calculation and will be presented elsewhere. It is a helpful 
feature that the zero syndrome, corresponding to no error, can be measured reliably. 

IV. SIDEBAND COOLING 

So far we have established that quantum information processing experiments can be carried out in the ground state 
of ions such as ^°Ca+, with good final readout. The choice of state avoids the need for very highly stabilised lasers, 
and the Zeeman shelving method avoids the need for very efficient fluorescence detection, making such experiments 
signiflcantly simpler than they would be otherwise. In this section we consider a further simplication, in which the 
need for very tight confinement in the ion trap is avoided. 

Cooling of the motional degrees of freedom is one of the most difficult parts of quantum processing experiments. So 
far only single ions have been cooled to the quantum ground state of a trap, first in one dimension only p9| ], then in 
three dimensions JT^ . Both experiments used sideband cooling P,pO| in the Lamb-Dicke limit. The Lamb-Dicke limit 
is the condition that the ion's motion is confined to a region small compared to the wavelength of the radiation under 
consideration. Since we are concerned with motional states near the ground state, of dimension ao, this condition is 
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rj = 27rao/A <C 1. From equation (|T|), this places a constraint on the confinement provided by the trap, parametrised 
by the vibrational frequency along the z axis, lOz S> Eji/h. Single ions of the species we are considering have recoil 
frequencies Eu/h in the region 6 to 30 kHz. Taking as an example a string of 10 ions we then require ~ 27r x 300 
kHz. To make the ions lie along a line, we ned also ujz ^ i^x.y ~ 27r x 3 MHz. This degree of confinement is 
experimentally accessible, but requires either very small closely spaced trap electrodes, or careful high voltage r.f. 
design. In the former case the coupling between the ion motion and noise voltages in the electrodes is enhanced, 
and in the latter the r.f. voltages are liable to be more noisy. Therefore there are advantages in using a less tightly 
confining trap. 

Another reason to avoid the Lamb-Dicke limit is that we need to drive vibrational-state-changing transitions (An — 
±1, where n is the vibrational quantum number of the fundamental vibrational mode) using the laser radiation, for 
information processing in the ions. A small value of 77 is not optimal for this purpose because it increases the off- 
resonant excitation of An = transitions when An = ±1 transition are driven. One could achieve the right conditions 
by first cooling to the ground state in the Lamb-Dicke limit, and then adiabatically opening the trap to 77 ^ 1 before 
processing, but it is interesting to know if this adiabatic opening can be avoided. Furthermore, for 77 ~ 1 the cooling 
rate is higher than for 77 ^ 1, which permits a lower steady-state temperature to be reached since heating effects 
(other than the random recoil from spontaneous emission) typically scale more slowly with trap tightness than the 
sideband cooling rate. 

Our approach to sideband cooling in the case of intermediate confinement, i.e. ~ 1, is essentially to use both 
the first and a higher order sideband. Second sideband cooling in a standing wave was considered by de Matos Filho 
and Vogel ||2l|, who showed that nonclassical motion can be obtained. We consider travelling wave radiation and are 
interested merely in obtaining low temperatures. Morigi et al. [ p2| recently discussed similar ideas to ours. They 
considered both r/ ~ 1 and the case 77 > 3 where new features become apparent. We provide a more complete analysis 
of the region < r/ < 2. Other theoretical analyses p3|-p7t have been restricted to first-sideband cooling and the 
Lamb-Dicke limit, though typically with a more thorough description of the atom-light interaction or the trapping 
potential or both. 

To describe the atom-laser interaction in the case of Raman transitions we adopt the standard theoretical device of 
considering the ground state manifold \M ~ ±1/2)^ as an effective two- level atomic system [ p8|j29| ]. This approxima- 
tion is valid when the excited state populations are negligible. The effective Rabi frequency describing the atom-laser 
coupling is n = f2_f2+/2A where ri± are the single-photon Rabi frequencies for excitation from \M = ±1/2)^ to 
an excited state such as P1/2 or Pz/2 (including Clebsch-Gordan coefficients), and A ^ il-t is the detuning from 
resonance of these transitions. The effective detuning is 

<5 = -.--c..+ -c.+_-M_dil±^, (14) 

where lol± are the laser frequencies and u)^ = jibB /h is the Zeeman splitting of the |M = ±1/2)^ levels. The effective 

linewidth F is equal to the optical pumping rate from = +1/2)^ to \M — —1/2)^ caused by a cr~ -polarized laser 
resonant with a transition such as \M = +1/2) ^ \Pi/2,M = —1/2). This laser radiation is either very weak and 
continuously present, or moderately strong and pulsed on when a spontaneous transition is desired. 



A. First sideband cooling 



We will discuss cooling of a single ion along one dimension. The generalisation to many ions and three dimensions 
is straightforward. Following Wineland and Itano |2^ , we will first account for the average change in total energy of 
the ion during one complete cycle of 'excitation' (|M = —1/2)^ \M = +1/2)^) followed by 'spontaneous emission' 
(|M = -1-1/2)^ — » \M = —1/2)^). The average is taken over many such cycles, in which the direction of spontaneously 
emitted photons varies randomly. The total state of the ion will be written \M — ±1/2)^ |n) where the ± indicates 
the internal state, and the integer n > indicates the vibrational state (external energy eigenstate) of the ion in the 
trap. 

Let 

/^„^|(/|e^''-|n)f . (15) 

where z is the position of the centre of mass of the ion. This is the relative strength of the different sideband 
components in the ion-radiation field interaction Using k • z = r]{a) + a) where a \n) — ^Jn \n — 1) one obtains 
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rain{n,f) Tl™ 2™ 

(/|e''"^|n) =e-'''/2y;^(i^)l/-"l y ^ ^. (i6) 

^ ' ' ' V J V /; to!(to+ 1/ - n|)!(min(n, /) - m)! ^ ^ 

The optical cross section for the 'absorption' transition |M = —1/2)^ |rt) \M = +1/2)^ |/) is 

afniS) - (To//«5 (-^ - [Ef - E,,)/h) (17) 

where g = gTTr/2 is the lineshape function (|^) normalised such that ^(0) = 1. If the occupation probability of 
vibrational state n is P„, then the rate of change in total energy of the ion during absorption, averaged over many 
cycles, is 

^ - E^"E - iI/^^+-)<^fn, (18) 

" / 

where / is the effective laser intensity. 

The average change in ion energy during spontaneous emission is 



£^sp = f dQ^P{k) J2 {En - Ef) I„f 

= f dQ^P{k) J2 (/I e''^-^ \n) {n\ [ff^'^*, e"*"^-^] |/) 



= Er. 

Note that since (/| p |/) — this answer is independent of the distribution P{k) of directions of spontaneous emission, 
which is counterintuitive. If, instead, the absorption and spontaneous emission are treated as a single scattering 
process (see eq. (p7|)), then the same result for the mean energy change is obtained as long as the distribution P(k) 
is symmetric about the origin, as it is in practice. 

Combining the contributions of absorption and spontaneous emission, the mean rate of change of total energy, per 
cycle of absorption followed by spontaneous emission, is 

^ = ^ E ^" E (^/ - ^" + {S - {Ef - E^)/h) (19) 

^ " / 

This equation was derived in ]2^ as part of a more general discussion. The above shows that its derivation can 
be simple and physically intuitive. It is valid to separate absorption and spontaneous emission because ty pically 
the experimental procedure is to switch on the Raman lasers and optical pumping laser at separate times p9| , [l4| . 
However, later we will not separate them but treat the emission and absorption as a photon scattering process. 

Although we are interested in the case 77 ~ 1, it is instructive first to examine the Lamb-Dicke limit 77 <C 1. This is 
done by a series expansion of eq. ( p^ ) in powers of 77^. The expansion of //„ (eq. (p^) to order ry'* is 

Ifn = 5f.n - ??'(2n + 1) + ^(3n2 + 3n + 1)^ 

+5f,n-irfn{l - rj^n) + 5f^n+irf{n + 1)(1 - ?7^(n + 1)) 

4 4 
+Sf,n-2^n{n - 1) + Sf^n+2^{n + l)(?^ + 2). 

This can be obtained either from eq. ( |l^ ) or by expanding the operator e^^'^ in eq. (|l5|). The terms in eq. ( p^ up 
to order 77^ are 



d{E) I<JoEjiy-^ r,_ ^ , 2/~ ~ \ 

~(tt~ ^ huj+ ^ ^^^'^ 9m+l) + V (5m+2 - 9 m) 



+n {{gm+l - gm-l) + Tf{grn-l - Hm+l + .gm-2/2 - 2gm + 35m+2/2)) 

+7i^7?^(-5„i_2/2 - (jrn+l + 9^-1 + gm+2/2)} , (20) 
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where m = 1, 2, . . . indicates the sideband chosen for cooHng, and gm = g{mu}z) = r^/(4m^a;^ + F^). 

For first sideband coohng, m = 1, in the limit of well-resolved sidebands we ignore terms of order F^/w^ compared 
to 1 to obtain, for the steady state distribution {d{E)/dt — 0), 

51 + .92 + ??'(ff3 - 5i) + (") (-1 + ri'') + (n2) r;^ = (21) 

In the Lamb-Dicke limit one finds that the rate of driving of population from one level to another is independent 
of the level number n, so the principle of detailed balance leads to the thermal distribution defined by P„+i/P„ = 
^-huj^/ksT ^ ^ populations are then 

P„ = (1 - (22) 
and the mean and mean square quantum number are 

/ > S ; 91 S(l + S) 

w = T^> (-')-(r^- (23) 

from which {n^') ~ 2 (n)^ + (n) . Substituting this in (|2^) and solving for n, ignoring terms of order rf or F^/cj^, we 
obtain 

, . gi+g2-7?^(gi-ff3) ^ (5 - 32r?V9) F^ 

~ 1 - 2772 ~ 16 (1 - 27^2) cj2 > 

This result is precise in the Lamb-Dicke limit with well-resolved sidebands. Figure 5 shows 1 — Pq as a function of 



772 as given by (24) (dashed line) and by a more exact treatment (full line) which will be described shortly. It is 
seen that ( p^ indicates the right qualitative behaviour for moderate values of 77, though it is not precise, owing to 
a slight departure from a thermal population distribution in the more exact treatment. The main conclusion is that 
first sideband cooling works very well for values of the Lamb-Dicke parameter up to 7/ ~ 0.6. It is noteworthy that 
the effective temperature given by fc^T = Yilij ^ j \vl(Pq j P\) is approximately equal to the recoil limit fc^T — En when 
77 = 0.5, F — Q.2ljJz, and is sub-recoil for smaller linewidths. 



B. Second sideband cooling 

For second sideband cooling, 771 = 2, eq. ( pO| ) can be usefully simplified by neglecting terms of order OI/rj^T"^ /loI). 
In other words, we make immediately the assumption of well-resolved sidebands. The steady state of eq. ( po|) is then 
given by 

{n") Y + - .93 - y ) - (52 + h) = 0. (25) 

There are two interesting regimes in which this result can be examined. First, for a very tight trap with 772 <c F2/ti;2^ 
the 772 terms are neglible and we have 



92 + 93 ^ 13 

gi - 33 ~ 32 ■ 



(26) 



This can be understood as the condition that all the population is driven to the ground and first excited vibrational 
levels, but that these are almost equally populated since the laser drives off-resonant transitions from each to the 
other at almost equal rates. 

Secondly, let us consider a trap that is less tight though the sidebands are still well resolved: F2/a;2 ^ 772 <c 1. 
Since the population in levels 7i > 1 is optically pumped efficiently towards n — 0,1, most of the population is in Pg 
and Pi, and therefore (^n'^) ~ Pi ~ (n). Using this to express (772^ in terms of {n) in eq. (p5|), the 7/2 terms cancel, 
so once again the mean quantum number is given by (^). This reasoning is valid as long as (77.^) — (77) <C 32 + ff3- 

In summary, cooling on the second sideband alone does not allow the ground state population to reach a value 
close to 1. Equation ( |2^ ) is nevertheless useful as a further test (in addition to (p^)) on the numerical treatment to 
be described in the next section. 
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C. Cooling on two sidebands 



We now turn our attention to 77 > 0.6. For single-frequency laser radiation on the mth sideband, the populations 
Pn can be deduced from a set of rate equations, in which the rate of transfer of population from level n to level / is 
given by 

2 



^fnim) = — I duN{u) 



^ LU^X^m- i + n) + iT/2 



(27) 



where N{u) is the angular distribution of spontaneous emission. Note that we are now treating the atom-laser 
interaction as a scattering process (still in the low intensity limit ^ F, cj^). If one wished to calculate the case that 
absorption and emission are separate processes (i.e. a driven transition followed by optical pumping), one would take 
the modulus square of each term in the sum over j, instead of the modulus square of the whole sum. We find however 
that the two cases give almost identical results, so we will restrict our discussion to the form (p7|). 

For dipole radiation in three dimensions one typically uses N{u) = |(1 + v?). To remain consistent with the 
one-dimensional approach we have adopted, we take instead N{u) = {5{u -I- 1) -f 5{u — l))/2, which corresponds 
physically to spontaneously emitted photons propagating along the z axis. We find that this correctly reproduces the 
Lamb-Dicke limit behaviour described in the previous sections. 

When 77 > 0.6, single frequency radiation will not cool the ion to the motional ground state. We therefore consider 
the next simplest case, namely two laser frequencies, one tuned to the first sideband, the other to the sideband m > 1 
(both are red sidebands, i.e. below the carrier frequency). To allow two laser frequencies in the population rate 



equations (27), we perform an incoherent sum, giving for the rate of transfer from n to f 

Tf,,^Tf„im)+aTf„{l). (28) 

This corresponds physically to the case that the radiation fields at the two frequencies are not simultaneously present, 
but pulsed on one after the other, and we consider the behaviour time-averaged over such switching. The parameter 
a is the ratio of the two laser intensities. 

The general physical picture of the behaviour is that the laser on the mth sideband pumps the population towards 
the lowest m vibrational levels, and cools overall, while the relatively low intensity first sideband laser has a small 
heating effect overall, but efficiently pumps the population from the lowest levels into the ground state. The rate of 
this pumping into the ground state is approximately a(l]^/F)/oo/oi = a{^^ /T)r]^e~^^ . 

In figures 6a and 6b we present the result of numerically extracting the steady-state populations, for a set of 
population rate equations using the rates (^). The figures show the value of 1 — Pq for 77^ in the range to 3, for 
two values of the linewidth, and for m = 1 to 4. The numerical calculations used a truncated set of levels, from the 
ground state n — to n — rimax- Even though the population is almost all in the low levels, it was found that a 
large number of levels (n^ax = 100) had to be included in order to calculate 1 — Pq. This is because for 77 > 0.6 the 
population distribution has a long oscillating 'tail' which falls off as a low power of n. As a result, the results in figs 
6a and 6b have only approximately 5% accuracy on the parts of the curves which rise steeply as a function of 

We found by numerical experiment that the ground state population does not depend strongly on the laser intensity 
ratio a, and that the optimal value is approximately given by a = 1/(377'^). This choice was adopted for all the results 
presented here. 

The main conclusions from figure 6 are as follows. We find that in the limit of well-resolved sidebands, 1 — Pq 
is roughly proportional to T'^/uj'^ for the whole range of rj values shown. The best sideband to choose, in order 
to maximise Pq at given 77, is higher than one might expect. A naive argument would suggest that since resonant 
absorption results in a kinetic energy decrease by mujz , and spontaneous emission results on average in a kinetic energy 
increase of Eji — rfuiz, the best choice is 777 equal to the smallest integer larger than rj^ , ignoring the small heating 
effect of the first sideband laser. In fact we find the optimum m is the smallest integer larger than approximately 
2772 + 0.5. 

The poor performance shown in fig. 6 around 77^ = 1.27 and t;^ = 2 is due to zeroes in namely hsiv^ = 
3 — VS) — and /i2(''7^ = 2) = 0. The first makes the depopulation of level 3 inefficient when ?Ti > 3, the second 
makes the depopulation of level 2 inefficient when 777, > 2. The structure in the m — 3 curve at 77^ ~ 1.5, 1.3, 1.2, ... is 
due to nearly coincident zeroes of /„_3,„ and In-i,m causing population to accumulate at 7i = 8, 9, 10, . . . respectively. 
Such coincidences only happen when both sidebands are odd, or both even, which is why similar structure does not 
appear in fig. 6 for ttt, = 2, 4. 

The method of cooling on a pair of sidebands will not work for large 77, because the first sideband laser cannot 
efficiently drive transitions from level 1 to 0. This was noted by Morigi et at We can roughly estimate the 

maximum 77 for which cooling to the ground state is possible as follows. 
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The equations are too complicated to solve analytically when rj is not small. In order to make rough esti- 
mates, we make several approximations. First approximate Tfn by ignoring the distribution N(u) and assuming the 
interference terms in the sum roughly cancel, leaving 

^f^-~^^ ^fj^jnign-j-m + affn-j-l). (29) 

We are interested in the case of low temperature, where most of the population is in the low-lying levels. In the limit 
that most of the population is contained in Pq and Pi, the ratio r = Pq/Pi can be estimated to be r ~ Foi/riQ. If r 
is large then Pq — 1, so cooling to the ground state has been achieved, and 

Po ^ . (30) 

i 01 + i 10 

We now examine the ratio 

Toi E^o lojIjiiOm+j-i + ag]) 



rio E^o hijoig-m+j + agj+i) 



(31) 



Since to > 1, the only resonant term in the two sums is the cjj term in the numerator for j = 0. Therefore the ratio 
may be written, in the limit <^ lj^, 

Toi /oo/oi + P _ e-^^\^+l3 

T^o ^ ^ " ~P ^ ^ 

where (3 is OiV"^ /ujD and we used the expression for /go and /qi given by (p^). The factor a cancels provided it is 
larger than ~ F^/w^. We estimate (3 as follows. By definition, 

oo 

- y^Jijijogj+i 

3=0 



^ j=0 ' 

where we used ( p^ ) for Iij and IjQ. The last step used — + which is true for large j and enables the 

sum to be performed. We find by numerical calculation that this overestimates /3 by a factor ^ rj^ . The sum in ( |3^ ) 
is equal to the Bessel function Io{2rj^) whose value for 2??^ > 2 is 

/o(2r;2) = ^=(l + 0((4,;)-2)). (34) 

Substituting this in (p3[) we obtain 



The most important feature of this calculation is that whereas /oo^oi falls exponentially with rf , we find that (3 does 
not, therefore it is impossible to achieve /oo^oi > P (and hence a large ground state population) once exceeds a value 
that depends insensitively on F/w^. In other words the optical pumping from n = 1 to n = becomes exponentially 
inefficient. We can estimate the limiting value of rj by setting /oo^oi = giving 

^ In {2uojT) + J In {if) + 2. (36) 

where the power •q~^ in equation ( |35| ) has been replaced by 77^^ in agreement with our numerical calculation. Equation 
(^) gives a reasonable estimate of the maximum value of rj that will permit accumulation of population in the ground 
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state by cooling on the first and a higher sideband. For F = O.luz, the hmiting value is 77 ~ 3, in agreement with the 
results of [^ . 

In conclusion, let us comment on the best choice of experimental parameters implied by this study. Our electron- 
shelving measurement measurement technique has as its major advantage the possibility of a high signal-to-noise 
ratio, since fluoresence can be driven for a long time before an unwanted transition makes the signal ambiguous. 
This renders information processing experiments in the ground state of '*''Ca'^ and similar ions feasible. Our study of 
sideband cooling has indicated the best choice of laser parameters in order to achieve a large population of the ground 
state of motion. In order to retain the simplicity of only cooling on a single sideband, combined with a not too great 
suppression of vibrational state-changing transitions, the choice rj = 0.5 is suitable. 
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FIGURE CAPTIONS 



Fig. 1. Low lying levels in ''Ca+. The Zeeman structure of the levels is shown, though the level spacings are not to 
scale. The wavelengths and inverse spontaneous rates of the electric-dipole allowed transitions are given. The circled 
numerals indicate the numbering of the levels adopted for equations (H) to (pT|). 




Fig. 2. Shelving efficiency e as a function of magnetic field in Tesla, for two values of the Rabi frequency of the 
shelving pulse: dashed line fiis = Ti^; full line Cli^ = ri5/4. The shelving pulse duration was chosen to maximise e 
for each value of magnetic field and Rabi frequency. 

Fig. 3. Shelving efficiency as a function of shelving pulse duration, for fli^ = Fis, B = 0.01 Tesla. Full line: 
solution of full rate equations; dashed line: equation (p^. 



Fig. 4. Steady state population of P1/2 level, N_y^ + N]^^^, during the detection phase, as a function of magnetic 
field. The observed fiuorescence signal is proportional to this population. The two laser beam intensities are set to 
give equal Rabi frequencies, f2 = = ^2a- Reading from the upper (full) curve to the lower (dash-dotted) curve, 
this Rabi frequency is r2 = 10r4, 5r4, 2r4, r4 respectively. 

Fig. 5. First sideband cooling. The figure shows the difference between the ground state population Po and 1, as 
a function of rf for first sideband cooling with F = O.lojj. Full line: numerical solution of population rate equations 
using rates r^„(l) (eq. (p^)); dashed line: analytic solution ( p4|) in Lamb-Dicke limit assuming a thermal population 
distribution, giving 1 — Pq = s ~ (n) . 

Fig. 6. Cooling on two sidebands. The curves show the difference between Pq and 1, as a function of 77^, for cooling 
on the first sideband and the mth, where m = 1 (full line), 2 (dashes), 3 (dots), and 4 (dash-dot), with a = 1/(37^^) 
which roughly maximises Pq. (a): F = O.lcj^; (b): F — Q.2lOz- The curves are obtained by numerical solution of the 
steady-state population rate equations with 100 vibrational levels included in the calculation. 
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